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Abstract. 

The Fokker-Planck equation for a heavy particle in a granular fluid is derived from 
the Liouville equation. The host fluid is assumed to be in its homogeneous cooling state 
and all interactions are idealized as smooth, inelastic hard spheres. The similarities 
and differences between the Fokker-Planck equation for elastic and inelastic collisions 
are discussed in detail. Although the fluctuation-dissipation relation is violated and 
the reference fluid is time dependent, it is shown that diffusion occurs at long times 
for a wide class of initial conditions. The results presented here generalize previous 
results based on the Boltzmann-Lorentz equation to higher densities. 
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1. Introduction 

In the hundred years since Einstein's remarkable paper on Brownian motion pQ, 
its impact has been felt on a wide range of problems in Physics, Chemistry, and 
Mathematics. The subsequent work of Smoluchowski [2j and Langevin [2j extended 
his physical model to more precise mathematical structures in the theory of stochastic 
processes. Attempts to recover these as idealizations of results from the actual 
Newtonian dynamics of the fluid came much later with the developments of non- 
equilibrium statistical mechanics [3] . The detailed description of Brownian motion as a 
function of particle size from first principles remains an open problem p. . A simpler and 
more controlled problem is that of a heavy atomic (small size) impurity in an equilibrium 
host fluid. Since there is a small parameter, the ratio of the host fluid particle mass to 
that of the impurity (and also small ratio of corresponding mass densities), a systematic 
analysis would appear to be straightforward. However, even in this simple 
perturbative analysis is found to be limited in accuracy to its leading asymptotic term 
and secular terms must be avoided at higher orders jSj. The asymptotic analysis for an 
impurity in a normal fluid at equilibrium leads to the familiar Fokker-Planck equation 
|H Ej. The objective here is to obtain the extension of that result for an impurity in 
a corresponding granular fluid. These results generalize earlier studies at low density 
Pi. 

The idealized prototypical model for a granular fluid is a system of smooth, inelastic 
hard spheres. This will be the model considered here. The inelasticity of collisions 
implies a continual loss of energy on each collision. Consequently, the usual equilibrium 
state for an isolated molecular system is replaced by a homogeneous "cooling" state 
(HCS), which serves as the environment in which the motion of a heavy impurity is 
considered. Although this reference HCS is inherently time dependent, an equivalent 
stationary representation can be given using appropriate velocity scaling. The time 
dependence of the probability density for the impurity particle in its phase space of 
position and velocity, can then be analyzed by methods similar to those for normal fluids. 
For example, the general case of diffusion of an impurity particle of arbitrary mass in 
a fluid undergoing HCS has been given recently in ^Uj. In this work, a Zwanzig-Mori 
projection operator method ^Tj is used (see | Appendix A ) to obtain a formally exact 



kinetic equation vaild for this general case. The Fokker-Planck equation then follows 
from this kinetic equation to leading order in the ratio of a fluid particle mass to the 
impurity mass. The analysis is complicated by the singular nature of the hard sphere 
interaction, but this mass ratio expansion has been discussed in detail for hard, elastic 
collisions ^2] ■ An interesting new complication for inelastic collisions is a restriction on 
the cooling rate of the host fluid relative to the collision frequency of the impurity. If 
this is too large, the expansion is no longer valid. This point has already been noted 
and discussed elsewhere [KJ US] • 

The rest of the presentation focuses on the similarities and differences between 
impurity motion in normal and granular fluids. First it is shown that there is an 



Brownian Motion in a Granular Fluid 



3 



exact mapping of the granular Fokker-Planck equation onto that for elastic collisions, 
in appropriate dimensionless variables. Thus the physical mechanisms of fast velocity 
relaxation to a Maxwellian for homogeneous initial states, and fast relaxation to a 
diffusive equation for inhomogeneous states occurs for granular gases as well. In 
particular, the transition to a hydrodynamic stage (diffusion) takes place for all degrees 
of inelasticity. However, in terms of physical variables there are significant differences. 
One of these is the violation of equipartition, since the kinetic temperature of the 
impurity is different from that of the fluid in the stationary state. Also, the usual 
fluctuation-dissipation relation between the drift and diffusion coefficients no longer 
holds, as expected for a non-equilibrium state of the host fluid. One consequence is 
that the mobility coefficient, measuring the response to an external force (for instance, 
electromagnetic or gravitational), is not simply related to the diffusion coefficient by the 
usual Einstein formula. 

The Liouville equation for a granular fluid [TJJ[T3] is introduced in the next Section. 
The special HCS for an isolated system is described and its time dependence is removed 
through the introduction of dimensionless variables. It is shown that the HCS is 
a stationary state for the Liouville equation in this representation. This stationary 
representation is used to describe the motion of an impurity in the fluid. A formally 
exact kinetic equation is derived for the impurity particle in Sec El and the condition for 
its stationary solution is described. Next, in Section 0] the form of this kinetic equation 
is simplified to that of the usual Fokker-Planck equation JS] in the asymptotic limit of 
small ratio of fluid mass to impurity mass. The details of the reduction are given in the 
Appendices. Some of the most important similarities and differences between impurity 
motion in normal and granular fluids are listed in Section and summary remarks are 
given in the last Section of the paper. 

2. The Homogeneous Cooling State and its Stationary Representation 

The basis for the analysis to be carried out in this paper is the Liouville equation for a 
fluid of N smooth, inelastic hard spheres (d = 3) or disks (d = 2) of mass m and diameter 
a, and one impurity of mass mo and diameter cq. The formal structure of this equation 
has been discussed in detail elsewhere in the general context of impurity diffusion [10J, 
and only an overview will be given here. The position and velocity coordinates of the 
N equal particles will be denoted by vf, i = 1, • • • , N}, while those of the impurity 
particle by qo,VQ. The expectation value for some observable A(T) defined over the 
N + 1 particle phase space T = {q , q 1: • ■ ■ , q N , v , Vi,---, v N }, is given by 



where p(T, t) is the probability density for the state of the system at time t and the 
dynamics is determined from the Liouville equation 




(1) 



(d t + L)p(T,t)=0. 



(2) 
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The generator for the inelastic hard particle dynamics is 

N 

L =L f + v - V -^T_((M), 

i=i 

N 1 ff JV 

L f = E^-^-jEE^j)- ( 3 ) 

?=i ?=i 

Here, L/ is the generator for the dynamics of the fluid particles alone and the subscripts 
or labels refer to the impurity particle. The second term in the expression of L 
generates the free streaming of the impurity. The operators T-(i,j) and T_(0, i) 
describe binary scattering between fluid particles and between the impurity and fluid 
particles, respectively, 

T_(i,j) = a d ~ l J dQQ( gij ■ ■ a\[a- 2 5(q i3 - <r)b£ - 5(q tJ + a)], (4) 

T_(0,i) = a d ~ 1 J dne(g 0l ■ a)\g 0t ■ a\[a^ 2 5(q 0l -W)b l - S(q 0i + (5) 

where dfl is the solid angle element for the unit vector a, cr = aa, is the Heaviside 
step function, and gjj = t>j — Vj and q^ = q« — q^ are the relative velocities and 
positions, respectively. Moreover, W = aa with a = (a + a ) /2 . The operators b^ 1 
and 6q/ determine the pre-collision velocities in a restituting collision. They are defined 
by 

1 + a, 

~2a 
1 + a 

and 



b tj Vi=Vi- -7^(gij ■ 
1+a 

K/ V i = v j + -^-(gij • ^i 5 "- ( 6 ) 



(l + a )A 
v o = v o [9oi ■ tr)a, 

,-1 (1 + a )(l - A) 

«o 

We have introduced A = mj (mo + m), that is the small parameter of the system 
in the Brownian limit to be discussed in the following. In addition, a and ao are the 
coefficients of normal restitution for fluid-fluid and impurity-fluid collisions, respectively. 
These coefficients measure the degree of inelasticity of collisions and take values in the 
interval < a, a < 1. 

Average kinetic (granular) temperatures for the fluid, T(t), and the impurity 
particle, T (t), are defined by 

1 N 1 
m^—Y^m^^-mv^t), (8) 

i=i 

T (t) = ^m (v 2 ;t) = i mo ^(t). (9) 
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The average velocities v(t) and v (t) introduced above will be called thermal velocities 
of the fluid and the impurity, respectively. Since all the collisions are inelastic and the 
effect of the impurity on the fluid is negligible in the limit of large N, T(t) decreases 
monotonically in time for an isolated system with a < 1. More precisely, 

d t T(t) = -C(t)T(t), (10) 

where ((t) is the "cooling" rate of the fluid due to inelastic collisions, 

C(*) = -^77w / dTvlT4l,2)p(T,t) > 0. (11) 



t?(t)d 

Similarly, the time evolution of T is given by 

d t T (t) = -Co(t)T (t) (12) 

with 

2N f — 

Co(*) = -=277^ / dTv 2 T40,l)p(T,t), (13) 



but it can not be concluded that the impurity temperature is monotonically decreasing, 
in general. Its behavior depends in detail on its initial value relative to that of the fluid 
particles. 

Since the fluid is cooling, there is no stationary solution to the Liouville equation 
for an isolated system. Instead, it is postulated that there is a solution whose time 
dependence occurs entirely through T(t) and T (t), having the scaling property 

Phcs {v,t) = [B 0>hcs (t)r d [B hcs (t)r Nd pIcs (H-X,q; X}) , (14) 

with the dimensionless variables defined by 

Q* = Qi/A Qo = <W^ v * = Vi/v hcs (t), Vq = v /vo,hcM- ( 15 ) 

The coordinates have been scaled relative to £ = (na^" 1 ) -1 , where n is the number 
density of particles, so that I is proportional to the mean free path. The subscripts 
hcs on the thermal velocities denote their values calculated for this particular solution. 
The dimensionless distribution function p* hcs is invariant under space translations and 
therefore Ph cs (T,t) represents a spatially homogeneous state of the fluid plus the 
impurity. The scaling form in equation ffT^j) implies certain constraints on the cooling 
rates. Let us introduce reduced cooling rates by 

c "- = iut)' c °'- [7(t)l= iw^r (16) 

From Eqs. (fTTj) and (fT¥j) it is easily seen that the reduced cooling rate of the fluid, Q cs , 
is time independent. Moreover, the cooling rate for the impurity in this state, Codes' 
depends on time only through the ratio of temperatures 



7(t) = ^T = Vn^- (1T) 

Similarly, it is found that the Liouville equation for p* hcs in these dimensionless variables 
depends explicitly on time only through j(t). Since p* hcs is time independent, solutions 
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exist only if ^(t) is a constant. The cooling rate equations (|1U|) and (|12j) then give 
immediately the requirement 

(L = iQms (7) • (18) 

In the particular state we are considering, the fluid particles and the impurity 
particle cool at the same rate, i.e. (hcs(t) = Co,/ics(^)- This should be viewed as a 
condition for the existence itself of this state as it fixes one of the two temperatures in 
terms of the other. However, this does not imply of course that the two temperatures 
are the same. In fact, it is found that they must be different except in the case of 
mechanically identical particles, as discussed below. With these results, the equation 
for p* hcs becomes independent of time and has the form: 

\&» E ^ • WPh~) + r (7) Pics = 0, (19) 
i=0 1 

where L* is the Liouville operator in dimensionless form (see below). The self-consistent 
solution to the coupled set of equations formed by equation (fTTIj) and the cooling 
equations determines the homogeneous cooling state (HCS). It is the analogue of the 
Gibbs state for molecular systems and reduces to it for a = a = 1. 

In dimensionless form, as defined by p* hcs1 the HCS is time independent. This 
suggests a transformation of the Liouville equation to dimensionless variables such that 
the HCS is a stationary solution in that representation. This can be formally carried 
out in the following way. We are going to scale the positions and velocities as given 
by Eqs. (fT5|) . where now Vh cs (t) is defined by Vh cs (t) = (2T/ lcs /m) 1 / 2 and Th cs {t) is the 
solution of the equation 

d t T hcs (t) = -Chcs(t)T hcs (t), (20) 

with the initial condition T hcs (t ) = T(t ), t being some arbitrary time. Of course, since 
in the general case the system is not in the HCS, it is T(t) 7^ T hcs {t). Moreover, in the 
scaling vo^csif) is chosen as given by vq^cs = IVhcsit), where 7 is the time-independent 
parameter identified by equation (fTHj). The appropriate dimensionless time scale s is 
proportional to the average accumulated number of number of collisions for the fluid 
particles in the reference HCS, 

s(t,t )= fdt'^fl. (21) 

For a general state whose distribution function is p(T,t), the dimensionless distribution 
p*(T*, s) is defined by 

p(T,t) = [B 0>hcs (t)}- d [B hcs (t)r Nd p*(T*,s), (22) 

where T* = {<7q, • • • , q* N) Vq, ■ ■ ■ , v* N }. Substitution of this into equation (j2J) gives the 
dimensionless Liouville equation 

(d s + C*)p*(T\s)=0, (23) 
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with the operator C* defined by 



N 



-r 

2 >>hcs 



= 1 
N 

E 







_d_ 
dv- 



(24) 



for arbitrary $(T*). The detailed form of the dimensionless binary collision operators 
T*(7;0,i) is given in | Appendix A The first three terms on the right side of equation 
(|24|) describe the dynamics of the impurity particle. Since the derivative with respect 
to s in (|2*H|) is taken at constant v$ and v*, some effects of cooling become explicit in 
the terms proportional to the cooling rate. The last two terms on the right-hand-side of 
()24j) generate the dynamics of the fluid without the impurity, again with explicit effects 
of cooling. As a consequence of these cooling terms, it is seen that the consistency 
condition for a stationary solution is the same as equation (|19|) for the HCS solution, 

£*plcs = 0, (25) 
since in equation (fiT?|) it is: 



N 



L*( 7 ) = L* f + 1V * ■ V - 7 E T *-(T, 0, i). 



i 



(26) 



Therefore, the dimensionless form of the Liouville equation, equation (|24j) . will be 
referred to as the stationary representation since it supports a stationary state and 
that state is the HCS. 



3. Kinetic Equation for the Impurity Particle 

Consider now the case of the fluid and the impurity in the HCS. At some instant (taken 
to be t = t = 0) the impurity is observed to be at position R with velocity Vo- This 
initial state can be represented as 

p (Y Q)= Pfecs (T) 5 (R - q ) 5 (V - v ) 

J alTp hcs (r) 5 (Rq - q ) 5 (V - v ) ' 

This is equivalent to 

* / r * n \ _ Phcs (£*) j (^0 ~ Qo) j (^0* ~ ^o) __ n * f-p* y*\ ( q\ 

The second equality introduces a convenient notation for the state in which the impurity 
variables X = {R, V} are sharply defined. In the remainder of this discussion, the 
stationary representation in the dimensionless units will be used, although for simplicity 
of notation the asterisks will be omitted. The probability F(R, V, s) for the impurity 
to have position and velocity R, V, at "time" s is then 

F(X, s) = j dT 8{X - x Q )p(T, t) = J dT5 (X - x ) e~ Cs p hcs (r ; X ) , (29) 
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where, trivially, Xq = {^o^o}- A formally exact kinetic equation for this probability is 
obtained in |Appendix A| It has the form 



ds 1 dR 1 2 dv 



F(X,s) 

+ J ds J dX'M(X, X\ s - s')F(X\ s) = 0. (30) 

The linear operator A, defined in equation (jA.16)l in |Appendix A[ describes the exact 
short time effect of collisions, including the initial correlations among particles in the 
HCS. If velocity correlations are neglected, A reduces to the Enskog-Lorentz collision 
operator. The form of the operator M(X, X', s) is given by equation ()A.27j) . It describes 
dynamically correlated collisions of the impurity with two or more particles that develop 
in time. These contributions are vanishingly small at short times and generally negligible 
at low densities. However, at very high densities, they describe the dominant effects of 
"caging" as the impurity particle becomes more localized due to collisions. 
The stationary HCS, Fh cs (V), is determined from equation (|3T)|) by 

Cj fW ' ^ VF ^ V ^ + 7AF to (V) = 0. (31) 
Use has been made of the property given by equation ()A.18|) . i.e. 

J dX'M(X, X\ s)F hcs {V) = 0. (32) 

Thus, the HCS is a stationary point of M but not of A. Further comment on this is 
given in the next section. Also, the particularization for the HCS of the expression for 
the cooling rate given by equation (fT3*j) can be expressed in terms of F hcs as 

2W f 

(o,hcs = — / dv v Q AF hcs (v ), (33) 

W being the volume of the system. The results of this Section are still exact and serve 
as a suitable starting point for the small mass ratio expansions to be discussed in the 
next Section. 

4. Fokker-Planck Limit 

In this Section, it will be shown that the kinetic equation (|3U|) reduces to a Fokker-Planck 
equation for asymptotically small m/m . The analysis proceeds in the following way. 
First, it is shown in |Appendix B| that the operator A becomes in the limit m/nio — > 

7 AF(X, s) = -Vo4? ■ (v + ^4t) F (X, s), (34) 



dV V IdV, 

where To and G are collision integrals, independent of time and velocity, given by Eqs. 
(|FHTJ) and flE20j), respectively, 
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G= 7 -*(l + ao )A C2| (36) 

C\ and C2 being real, positive numbers, of the order of unity. Note that in the limit we 
are considering, 7 _2 A is a constant, independent of the mass ratio. For elastic collisions, 
G — > 1 and A becomes the usual Fokker-Planck operator in dimensionless form. The fact 
that G 7^ 1 here can be considered as a violation of a fluctuation-dissipation relation. 

4-1. Impurity HCS 

Let us analyze the consequences of this for the impurity HCS. Use of equation (}34j) into 
equation (J3T]) gives immediately 



d 

dV 



e)V + lG- 9 



2 dV 



hcs 



0, (37) 



where we have introduced the parameter 

e = ^. (38) 

It is seen that the solution to this equation is a Gaussian. Normalization and the 
definition of the impurity temperature require that 

F hcs (V) = TT- d ' 2 e- v2 . (39) 

Then equation (|3*7)l gives the relation 

G=l-e. (40) 

This allows determination of the impurity temperature in the HCS in terms of the fluid 
temperature. The latter is determined from equation (fTUj) using for ( hcs the expression 
identified bellow, and has a universal value, depending only on a and being independent 
of the initial conditions at long times. The impurity temperature in the HCS is then 
determined from equation (jlUj) with G given by equation (JSHJ, 

T„(«) = (1 0^T(t), (41) 

following that T (t) ^ T(t). 

A further consequence of the above relationship is the condition e < 1, i.e. 

^ < 1. (42) 

Thus the small mass ratio limit for granular fluids entails, in addition, a restriction on 
the inelasticity of the fluid-fluid particle collisions. This follows because (hcs oc 1 — a 
and r oc A, so the fluid inelasticity must decrease as the mass ratio goes to zero. 
Asymptotically, for small 1 — a it is found that (see | Appendix C ) 



4 vr (^-i)/2 

a -~viTO (1 - Q)xW ' (43) 
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where x( a ) is the fluid-fluid pair correlation function for particles at contact and T(y) 
is the Euler Gamma function. Then the condition (14*21) becomes 



A^)/\{a){\ - a) ^\ rf - 1 <1 (44) 



V2T(d/2)Cid(l + a )A \W, 
In practice, for any choice of A << 1 the value of a must verify 1 — a << 1, moreover 
being consistent with equation (J44|) . Otherwise, the derivation of the Fokker-Planck 
limit for A does not hold. The correct limit for e > 1 is described in reference |13j . 
Equation (J44j) can be made more explicit in the low density limit of the fluid. Then, 
the correlation function appearing in the expression of Ci, equation (jB.18|) . can be 

(2) 

set equal to unity, and g\ can be neglected. In this way, it is obtained 



a fa 



-i 



. . < 1, (45) 
v^(l + a )A \WJ 

that agrees with the result derived from the Lorentz-Boltzmann equation reported in 
reference jH]- 

In summary, the stationary state condition fixes the form of A as a Fokker-Planck 
operator with friction constant Tq and diffusion coefficient (1 — e), 



AF = -7-^0 ° 



dV 



V + ift *° 



(46) 



2 V ' dV 

where e has been defined in equation (|3*fi|h Moreover, the stationary state equation (pTfj) 
simplifies to 

(1 - e) J- - (v + 14t) F UV) = 0. (47) 



dV V 29V, 
4-2. Fokker-Planck Equation 

With the form for A given in equation (|46p. the kinetic equation ()30p becomes 
(d s + 7 v • V R ) F(X, s) = (1 - e) r A • (v + i F(X, s) 

+ y ds' y dX'M(X, X\ s - s')F(X', s'). (48) 

The first term on the right side is now a Fokker-Planck operator of the usual form, but 
the "friction constant" r obtained from A, has been renormalized by (1 — e), due to 
the additional operator representing cooling. 

The values of 7, e, and r are now known and it remains only to determine the 
asymptotic form of the last term on the right hand side of equation ([48)1 . There are two 
steps in that process. First, since M(X, X', s') involves two binary collision operators it 
turns out to be proportional to (A 7 -1 ) 2 . By exploiting this, in |Appendix B it is found 
that 

y dX'M(X,X',s-s')F(X,s') = A. (v + lJL^g c{s - s > )F ( X ,s'), (49) 
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where Q c (s) is the correlation function 

g c {s) = 2 -(F + {s)-F„) } . (50) 

The brackets denote an equilibrium ensemble average for the fluid in the presence of 
an infinitely heavy impurity at rest. The phase functions F± are the total rate of 
momentum transfer by the fluid particles to the impurity, 

- r 

F± = 7" 1 (1 + af) Aa d - 1 ^ / dQQ(^ Vi ■ £)<5(q 0i + W)( Vi ■ BfB. (51) 

i=\ 

For continuous potentials they would be simply the total force on the impurity, but 
are determined here for hard spheres by the binary collision operators. The final step in 
the reduction of term involving M(X,X', s — s') is to recognize that the time scale for 
the impurity as described through F(X, s) is slower by a factor 7 than the decay time 
for the fluid correlation function Q c (s). To leading order in the mass ratio this gives 

f ds'g c (s-s')F(X\s')~T c F(X',s), (52) 
Jo 

with 

POO 

ds'GcW). (53) 



With all these results, the final form for the Fokker-Planck equation is obtained: 
(d s + 7 V ■ V R ) F(X, s) = (T - er ) A . (y + I F ( X , s), (54) 

where 

r = r + r c . (55) 

is the total friction constant. In this dimensionless form, the inelasticity of the granular 
fluid is suppressed and occurs only through the effective friction constant T — eT . 



5. Comparison of Granular and Normal Fluids 

The Fokker-Planck equation plays an important role in demonstrating explicitly some 
important general beliefs of non-equilibrium statistical mechanics. Among these is a two 
stage approach to equilibrium. In the present case of a granular fluid, the role of the 
equilibrium state is played by the HCS. The first stage is a fast relaxation of the velocity 
distribution on the time scale of a few collisions, followed by a slower transition to spatial 
homogeneity and full equilibrium by hydrodynamic processes. This follows directly from 
the Fokker-Planck equation, where the only hydrodynamic process is diffusion. 

The recovery of the usual Fokker-Planck equation in dimensionless variables, 
equation ([54)1. means that many of these qualitative concepts of Brownian motion for 
normal fluids can be transferred to granular fluids as well. The general solution (Green's 
function) for the Fokker-Planck equation can be constructed to demonstrate the two 
stage relaxation. This has been discussed in some detail for an analysis based on the 
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low density Boltzmann-Lorentz equation [10, and the analysis there applies here as well. 
The main similarities between normal and granular fluids are: 

• There is a universal stationary state, Gaussian in the velocity and spatially uniform, 
approached at long times for all physically relevant initial conditions. For initial 
homogeneous states, the approach to this stationary state is exponentially fast in 
the time scale s determined by the collision number. 

• The spectrum of the Fokker-Planck equation includes an isolated point representing 
diffusion. This is the smallest point in the spectrum, corresponding to the slowest 
possible excitation. The next slowest mode decays exponentially fast in s relative to 
the diffusion mode. Therefore, for general initial conditions, the diffusion equation 
dominates for sufficiently long times. The existence of diffusion and its dominance 
applies for all values of the restitution coefficient a . 

This result is relevant because the validity of a hydrodynamic description for 
granular fluids is not self-evident. Furthermore, if that validity is granted based on 
empirical or phenomenological grounds, it is often assumed to be limited to weak 
inelasticity. The Fokker-Planck limit considered here provides an unambiguous example 
of the existence of hydrodynamics even for strong inelasticity, and a description of the 
initial transient period leading up to its dominance. 

There are also significant differences between normal and granular fluids that 
have been suppressed by the use of dimensionless variables. Returning to the original 
variables, equation (j5^| takes the form 

(d t + V • V R ) F(X, t) = A • jr [T(t)\ Y + \v [T(t)\ A I F (X, t).(56) 

The first term of the brackets is known as the "drift" term and describes the deterministic 
dynamics of the particle. The friction coefficient T{t) is related with the dimensionless 
one T* appearing in equation (pjj) by T(t) = v(t)T*/£. The second term represents 
diffusion in velocity space due to the finite temperature of the host fluid and gives rise 
to all statistical properties of the dynamics. The velocity diffusion coefficient is given 
by 

V [T(t)\ = ^A {r [T(t)} - eT [T(t)]} . (57) 

For Brownian motion in a normal fluid, the corresponding relationship is T> = 
(2T/mo)r. This is known as a fluctuation-dissipation relation (velocity diffusion 
representing fluctuations, the friction constant T representing dissipation). This relation 
is violated for granular fluids in two ways. First, there is the additional term proportional 
to e due to the cooling of the host fluid. Second, the impurity particle temperature T (t) 
differs from that of the host fluid, as shown in equation (pITj) for the impurity HCS. A 
related consequence of cooling and the two different temperatures, is a violation of the 
Einstein relation between spatial diffusion and mobility. The latter can be determined 
from the Fokker-Planck equation by adding an applied external force and calculating 
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the response of the average velocity of the impurity particle. The resulting mobility 
coefficient has no simple relationship to the coefficient of diffusion determined from the 
mean square displacement of the particle. 

An obvious effect of the presence of a time dependent temperature in is a 

nonlinear change in the relevant time scale. This follows from the fact that the fluid- 
fluid collision frequency sets the physical time scale and this is proportional to y/T(t). 
Consequently, the dimensionless collision frequency s depends logarithmically on t. This 
renormalization of the time is important for identifying the various excitations. For 
example, the mean square displacement of the impurity becomes linear in the collision 
number time scale s, not in t. 

6. Discussion 

The Fokker-Planck equation derived here has several qualifications related to the small 
parameter m/m Q . In fact, there are other parameters of the problem that can become 
large such that their product with the small mass ratio compromises the validity of the 
expansion. If the size of the impurity becomes large relative to the host fluid particle 
size, the mass densities of the impurity and fluid particles can be come comparable. In 
fact, this is a common experimental condition. In this case, it has been suggested [3] 
that the Fokker-Planck equation must be modified to a non-Markovian form, since the 
time scale separation in the analysis of the correlated collisions operator M is no longer 
justified. 

A second possibility is that the ratio of the impurity temperature to the fluid 
particle temperature, T /T, becomes large. This occurs when the cooling rate for the 
fluid is larger than the effective impurity-fluid collision rate (violation of the condition 
(142)1 above). Since the latter decreases with the mass ratio, this requires that the host 
fluid particle must be less inelastic. When this is not the case, the derivation given here 
is not valid and the Fokker-Planck description does not apply. Instead, the impurity 
particle executes "ballistic" rather than diffusive motion [T3]. 

The hard sphere interactions lead to some important differences from the results for 
continuous potentials. The collision operator in equation (130)1 has both an instantaneous 
contribution, A, and one representing finite time correlated collisions, M. The first 
is analogous to the Boltzmann-Enskog-Lorentz collision operator. Accordingly, there 
are separate contributions from each of them to the friction coefficient, T and r c , 
respectively. This is puzzling since for continuous potentials the friction coefficient is 
given by a single Green-Kubo expression, in terms of the time integral of the force 
autocorrelation function. However, it has been shown recently that this autocorrelation 
function develops a singularity for steeply repulsive potentials, leading precisely to the 
instantaneous contribution T for hard spheres [T7j. The residual non-singular part is 
given by the "force autocorrelation function" (150)1 . This is a general feature of the hard 
sphere interaction and is not related to the inelasticity of the collision. 

It is interesting that the operator M describing correlated many-body collisions 
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plays no role in determining the HCS distribution for the impurity particle or its cooling 
rate in the HCS. On the other hand, it does have an important effect for the dynamics 
of deviations from the HCS, through the correlated part of the friction constant r c . 
However, the final temperature of the HCS is independent of this constant. 

Within the context described above, the results are exact and apply for both dense 
and dilute fluids. Explicit evaluation of the friction coefficient as a function of the 
density and restitution coefficients is a difficult many-body problem. The contribution 
To is similar to that from the Enskog kinetic theory and in fact reduces to the latter 
for elastic collisions. Here, it is necessary to understand the pair velocity correlations 
in the HCS before further simplification is possible. The contribution from correlated 
binary collisions T c is more complicated. Even at low density and elastic collisions, an 
infinite sequence of "ring" and "repeated ring" recollisions between fluid and impurity 
particles must be calculated when the size of the impurity becomes comparable to the 
mean free path [T5] . 

Certainly, the motion of an impurity in a granular fluid provides conceptual, 
computational, and experimental challenges for the second century after Einstein's 
initial insight. 
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Appendix A. Formal Kinetic Equation 

The projection operator method for normal fluids provides a means to write formally 
exact equations for observables of interest in terms of collision kernels. These kernels 
are then the appropriate objects to study by means of chosen approximation schemes 
[TT] . In this Appendix, this method is used to obtain the exact kinetic equation for the 
probability distribution of the impurity particle, defined in equation (|29p. The method 
has been discussed in detail in many places for normal fluids and extends directly to the 
granular fluid when working in the stationary representation. Only a short summary is 
given here. 

The projection operator over the "relevant" part is defined by 



for an arbitrary phase function $(r). The distribution function ph cs (T; X) is defined in 
equation (j2HJ)- It is easily verified that V actually has the projection property V 2 = V . 
Moreover, it is 




(A.l) 




(A.2) 
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Consequently, Pp(T,t) retains the relevant information in order to evaluate the 
probability distribution for the impurity. To obtain an equation for Pp(T,t), we 
start from the Liouville equation (j2Bj) • where the binary collision operator for collisions 
between the impurity and the fluid particles is 

T_( 7 ; 0, i) = a*- 1 Jdne ^ 0i ( 7 ) ■ a] |g w ( 7 ) ■ B\ [a 2 5(q 0i - W)b - - 5(q 0t . + if)] , (A.3) 
where 

9oi(l) = v o ~ 7~V (A.4) 
and the operator b 0i is defined by 

= [0oA7) ' o- o-, (A.5) 

-: — i (1 + a )(l - A) 7 r _ ^ ^ 

6oi «i = «i + - [flf«(7) • o-] o-. (A.6) 

We operate on both sides of the Liouville equation with the operator P and 
rearrange the result to get 

(d s + PCP) Pp (r, s) = -PCQp (V; s) , (A.7) 

with Q = 1 — P. Repeat this analysis by operating on the Liouville equation with Q 
and rearranging in a similar way 

(<9 S + QCQ) Qp (r, s) = - QZPp (r, s) . ( A.8) 

Equations (|A.7|) and (|A.8|) are a coupled pair of equations for Pp and Qp. They are easily 
solved with the initial conditions Pp(T;0) = p (T, 0) = p hcs (T;X ) and Qp(T;0) = 0, 
following from equation (|28|). with the result 

(d s + PCP) Pp (r, s) - f ds' P£e- {s ~ s ' )QZQ Q£Pp (T, s') = 0. (A.9) 

Jo 

Finally, making some of the projection operators explicit and using equation (|A.2|) . gives 
the desired kinetic equation 

d s F(X,s) + J dX'B(X, X')F{X',s) 

+ J ds' J dX'M(X, X',s- s')F(X', s') = 0, (A. 10) 

with 

B(X,X') = j dT5(X-x )Cp hcs (T;X') } (A.ll) 

M(X, X', s) = - J dT5 (X-x ) Ce- sQXQ QCp hcs (V; X') . (A.12) 

The kernel B(X,X') represents the instantaneous mean field effects of the fluid 
on the particle. This can be made more transparent by rewriting its expression in the 
following way: 

B{X,X') = J dTS(X-x )Z8(X'-x )p hcs (T)F^(V'), (A.13) 
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where 

F hcs (V) = JdT 5(X - x )p hcs (T) (A.14) 



is the impurity distribution in the HCS. Then, using equation ([24)1 . it is found 

5(X, X') = 7 | ( V ■ Vh + A) 5 (X'-X) + &&lJL.[V6 (X-X')\ } , (A. 15) 

with the linear operator A defined by 

AS {X'-X) = - J dx 1 T4 r ,X 1 x 1 )5(X'-X)FW(X,x 1 )F^(V). (A.16) 

Here 

F®(X, Xi ) = nJ dx J dx 2 ■ ■ ■ J dx N 5(x - X) Phcs (r) (A.17) 

is the reduced two particle distribution for the impurity and one fluid particle, and the 
operator T_(7; X, X\) is still defined by equation (|A.2|i . but replacing x by X. Equation 
(|30|) follows by substituting equation ()A.15|) into equation ([A.lOp . As discussed in the 
main text, A is closely related to the linearized Boltzmann-Enskog-Lorentz operator for 
a granular fluid in dimensionless form. 

The kernel M(X, X', s) is somewhat more complex. If the complete system is in 
the HCS, i.e. p(T, s) = Ph cs (T) and, consequently, F(X,s) = F hcs (V), it is 

J dX'M(X, X', s)F hcs (V) = 0, (A.18) 
for all s. This is a consequence of the property 

J dX'£{T) Phcs {T,X')F hcs {V') = C Phcs (T) = 0. (A.19) 
Simplifications of M occur by writing 

N 

C = C + C f - 1 J2 T -^', 0, i), (A.20) 

i=i 

where Co and Cf are the generators for the dynamics of the isolated impurity and the 



isolated fluid respectively. Their expressions are directly identified from equation 
It is easily verified that 

V(£ + £ f )Q = 0, (A.21) 

and equation (jA.12|) thus becomes 

M(X,X', S )= 7 / dT6(X-x )J2 T ~(r,0,i) 
J i=i 

xe- sQZQ Q£p hcs (T) 6 (X'-x ) F£(V). (A.22) 
Next, use the identity 

c Phcs (r) $(r) = [£ Phcs (r)] $(r) + Phcs (r) £_$(r) = Phcs (r) £_$(r), (A.23) 
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valid for arbitrary $(T). The new generator £_ is defined by 

£_$(r) 



N Chcs d 



i=l 

AT , JV JV AT 



i=l i=l j^i i=l 



$(r), (A.24) 



i=l i=l j^i 

with the modified binary collision operators 

T_(i, j) = a^ 1 /" dfi^qy - < T)6(g, J ■ ^gy • S| (fer.i _ 1) , (A.25) 

T_( 7 ; 0, •) = ^ J cm 5(q 0l - <T)0(g Ol • £)|g 0i ■ <x| (V - l) . (A.26) 

For systems with elastic collisions, £_ is the generator for the time reversed dynamics, 
and it plays a similar role for the HCS time correlation functions of granular systems. 
The detailed proof of equation (|A.24|) will be given elsewhere. With this identity, 
equation (jA.22|) can be rewritten as 

M(X,X',s) = - 7 2 / dT5 (X-x )J2T-(r,0,t)e- sQCQ Q Phcs (T) 
J i=i 

N 

xJ2 T -(rAj)S(X'-x )F^(V). (A.27) 
Appendix B. Mass Ratio Expansion 

The expressions for A and M are given in Eqs. (jA.16|) and ()A.27|) . respectively. With 
the appropriate choice of H(xq,xi)) both have the form 



Z = J dx 5(X-x ) J dx 1 T^(r,0,l)H(x o ,x 1 ), (B.l) 

that using the expression of the binary collision operator can be transformed into 

Z = Jdx J dx l 5{X-x )a d - 1 J d06(g Q1 • £)|g 01 - a\ 

x [a 2 5(q i - &)b ± ~ 5 (Qoi + &)] H(x , Xi) 

= J dx J dx 1 H(x ,x 1 )S (R-r )a d_1 J dtt 0(g ol ■ <r)|g 01 ■ cr|5(q i + a) 

x [5(V-b 01 v ) -5(V-v )] , (B.2) 

where froi^o is obtained from equation (|A.5|) . 

b iv = v - (1 + a ) A(g 01 • ct)^. (B.3) 

In the limit of a very massive impurity, it is A 7 _1 << 1, and S(V — froi^o) i n equation 
!2"|) can be expanded to second order yielding 

dxo / dxiH(xo, x\)5 (R — r ) 
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d 1 d d 

-A(x , xi) • —5 {V - v ) + -B(x , xi) : ^—^—5 {V - w ) 



x / dxi 
with the definitions 



dx 5(X — x ) 
d 



A(x , Xi) + — — — : B(x , arij 



i7(x ,xi), 



(B.4) 
(B.5) 



(9w 2 (9i; (9^o 

Ai(x ,xi) = (l + a )Aa d - 1 J ^0[g ol (7)-CT]5(q O i+^)[9oi(7)-^] 2 ^ 

^■(xo^O = [(1 + a )A] 2 W d - 1 J dO9[g 01 ( 7 ) • £]%oi + <r) |^oi(7) • ^l 3 ^^' • (B.6) 

Comparison of Eqs. (jB.l|) and (jB.4|) gives the relation 

d 



dx 1 T_('y;0,l)H(x o ,x 1 ) 



1 a 



dv 



dx\ A(x , xi)H(xq 1 x\) 



+ 2dv' I dx ^{ x ^ x i)H{xQ,x x ) 



(B.7) 



It should be noted that this is not the complete expansion in the mass ratio, since there 
remains a dependence on 7. 

Appendix B.l. Evaluation of A 

For the particular case of A, equation (jB.7|) gives 

AF(X,s) = - [ dx 1 T_( r ,X,x 1 )Fl l 2 J s (X,x 1 )F^ s (V)F(X,s) 



d 

dV 



A(X) + iA. B (X) 



2dV 



F(X,s) 



(B.8) 



with 



dx S(x -X) I dx l A i {x ,x 1 )Fl l 2{X,x 1 )F^{V) 



l + a )A^^ J dx 5(x -X) J dvnp hcs (vi) J dQQ(g 01 -a) 



x(9oi ■ $) 2 g {2) ((T,V,v 1 )d i 



(B.9) 



B i3 {X) 



,/,■„,*(,■„- A') I dx 1 B ij (x Q ,x 1 )F^(X,x 1 )F^(V) 



[(l + «o)A]M- 



a 



dx 5(x -X) J dvuphcsivt) J dttO(g 01 



x(g 01 -afg (2 \w,V,v 1 )d t a J . (B.10) 

In the above expressions, we have introduced iphcs(vi) = n" 1 fhcs(vi), where fhcs(vi) 
is the one-particle distribution of the fluid in the HCS and the correlation function 
g^f^x, V, Vi) between the impurity and the fluid particles, defined through 

F™(W,V,v x ) = f hcs (vi)F hcs (V)g^(W,V, Vl ). (B.ll) 
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Because of the spatial homogeneity and isotropy of the HCS, the correlation function 
must be a function of the scalars a, V, V\, V ■ tx, v\ ■ er, and V ■ v%, where the hats indicate 
unit vectors. Furthermore, the characteristic impurity velocity is a factor 7 smaller than 
the characteristic fluid velocity. Therefore, the correlation function may be expanded 
to first order as 

gW(W, V, Vl ) ~ gf\a, Vl ,v l -a)+jV- gf\a, v^v.-a). (B.12) 
Now equation (jB.9|) can be written 

A(x) = 1 - 1 r v, r = rf + rf, (B.13) 

where 



rf~(l + a )A 



a 



a 



d-l 



dv X (p hcs {vi) I dflQ(-v 1 ■ S)gf ] (a,v 1 ,v 1 ■ er)(V • a 



x 



7 - 1 y- 1 (vi • sy - 2(«i ■ $)(v ■ a) 



(B.14) 



The term of order 7 _1 is of odd parity in v x (change V\ — > —Vi, <r — > —a) and, therefore, 
it vanishes when integrated over the isotropic distribution ifihcs{vi). This leaves the 
dominant term as 

d-l 



-2(l+Qto)A 



dv^hcsivi) / dQQ(-v 1 -a) 



x(vi • <r)(V • <r) 2 g$ i) (<r,v 1 ,v 1 ■ a). 

(2) 

Similarly, the contribution to T from J to leading order is 
r<i ] = (1 + « ) A f^j J dvwncivi) 

Combining these gives 



(B.15) 



1 + ot ) A (a 



a 



d-l 



Ci 



m (1 + a ) fa 
tuq 2 V a 



d-l 



Ci, 



(B.16) 



(B.17) 



x 



(B.18) 



where C\ is a pure number, 

Ci = 4 / dw^hc^i) y tZne(-Vi ■ ■ ?)(V • a) 

'y ■ ^)9o\^,v 1 ,v l ■ a) + («i • CT)\>-g5 2) (CT,t;i,Vi ■ ct) 

The analysis of B^X) follows in the same way, with the dominant contribution 
being: 

B ij (X)=j- 1 Gr 5 ij , (B.19) 

where 

(B.20) 



G = 7' 2 (1 + go) A ^ 
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and C2 is another pure number, 

C * = ~C~~d J dvi ^ hcs ^ Vl "> J d^®^- a)(vi -afg^ (a, vi,vi- a), (B.21) 
so that G is independent of the mass ratio in the Brownian limit we are considering. 

Appendix B. 2. Evaluation of M 

To evaluate M(X,X',s) in the small mass ratio limit, first equation (|B.7j) is used into 
equation (|A.27J) to get 

f 8 N 

M(X, X', s) ~ - 7 2 / dT6 (X-x Q ) ^- ■ ^ A(x , xfc-^Qp^ (r) 

J i=i 

JV 

xJ2t-(t,o,j)Hx'-x )f- c ](v') 

g 

7 av F 

AT 

xJ2T-(r,0,3)HX'-xo)F^(V'), (B.22) 
where we have defined 

N N r 



J dT6 (X-x )F + (T)e~ sQCQ Qp hcs (T) 



F + (T)= 1 J2A(x ,x t )= 1 (l + a )Aa d - 1 J2 [ dne^ ■ a)5(q oi + W)(g c 

8=1 1=1 ^ 

~ 7 - 1 (l+a )A^- 1 J] y dfieC-Wi-^qoi + ff)^.^ 2 ?. (B.23) 



Upon writing equation (jB.22)) . we have neglected the term proportional to B in equation 
flB.7|) . because it is of the order 7 _3 A 2 , while A is of the order 7 _2 A. In the case of 
A, the situation was different because the leading contribution from A happened to 
identically vanish. 

The action of the operator T_(7;0, j) on 5(X'—xq) can be calculated to leading 
order for a similar analysis to the one carried out at the beginning of this Section, 
obtaining 

T_( 7 ;0,j)5(X'-x ) =^ d ~ 1 J dQ6(c loj -W)e(g or a)\g or a\ 
x 6(x'- fco/xo) - S (X' - x ) 

- - (1 + a ° )A ^' 1 / dns( qQj - w)e(g oj ■ ?)(g 0j • af 

«o J 

xa~6(X'-x ). (B.24) 

Then, 

M(X,X\a)~—!L. J dY5 (X - x ) F + (T)e- sQZQ Q Phcs (T) F_(T) 
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xF^(V')~6(X'-x ), (B.25) 

with 

N 



^ 7 (l + a )A _ d _ 1 yy /" doe , . .^(qoi+wKvi-a) 2 *. (B.26) 

Since both -F+(r) and -F_(r) are proportional to 7~ 1 A, M(X,X',s) is of order 
(7 _1 A) 2 ~ m/mQ. This is the same order in the mass ratio as the leading contribution 
to A. Therefore, all additional corrections on the mass ratio in equation (jB.25|) should 
be neglected by consistency. In particular this means we can substitute 

e -sQCQ Q = e -sQC Q _ e sQC' f Q = e -sc' fQ ( R27 ) 

where C'j is the generator for the fluid dynamics, including the interactions of the 
particles with the impurity at rest. More explicitly, 

N 

£/ = £/-XX-(M, (B-28) 

1=1 

where Cf is the generator for the isolated fluid and 

T'_(0,t) = a*' 1 J dne(-V! ■ a)\v! ■ a\ [a 2 5(q 0l - W)b'r l + 5(q 0l +W)} , (B.29) 

b'r'v, = Vi - ( Vi ■ (B.30) 

Since C does not contain dynamics for the impurity, equation (|B.25|) can be transformed 
into 

M(X,X',s) = ~- J dTF + (T)e- sZ 'f Phcs ^)F-^) 

x5(X - rro)F^(V) • ^5(X - X'), (B.31) 
where the relation Vphcs(j)F^(T)S(X' — x ) = has been used. Then, 

J dX'M(X, X\ s - s')F(X>, s') = A . Gc ( s - s ').(v + ± F(X, s'), (B.32) 
with 

G c (s) = 2 J dTF + {T)e- st f Phcs {T)F_{T)8{X-x Q )F^ s {V)- (B.33) 
Because of the homogeneity and isotropy of the HCS, it is 

G c (s) = g c (s)\, (B.34) 
where I is the unit tensor, and 

g c ( s ) = - [ dv f \f + {t) ■ e- sZ 'f Phcs (r f \x)F_(r)} . (B.35) 

a J 1 J x =x 
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Here V f is the phase space associated to the fluid particles and 

p(T f \X) = [ Phcs (r)L =x KcliV) (B.36) 

is the conditional HCS probability distribution of the fluid particles given that the phase 
space coordinates of the impurity are X = {R, V}. To leading order in the mass ratio, 
we can neglect the velocity correlations between the fluid particles and the impurity, 
as already discussed in the evaluation of A. Then, Q c {s) becomes independent of X. 
Moreover, in the limit a — > 1, ph cs {T f\X) can be substituted by its equilibrium form. 
In this way, equation ()B.32j) becomes 

J dX'M(X, X', s - s')F(X', s ') = ^-{v+ l -^jg c {s- s')F(X, s') . (B.37) 
This is the expression (J4"U|) quoted in the main text. 



Appendix C. Evaluation of the Cooling Rate of the fluid 

The dimensionless cooling rate for fluid in the HCS defined by Eqs. (fTU|) and (|TB|) can be 
reduced to an integral over the appropriate reduced distribution functions. The details 
are given in the Appendix A of reference JU]> with the result 

Chcs = ^ (1 - a 2 ) a d - 1 J d Vl J dv 2 J dQ /(2>( Vl , v 2 , q 12 = a) 

x6(-g 12 -£)|g 12 -£| 3 , (C.I) 
where the fluid two-particle distribution function, 

f^s(x 1 ,x 2 ) = N(N -l)Jdx Jdx 3 ---J dx NPhcs (T) , (C.2) 

has been introduced. Due to the condition e < 1, equation (f^)). the restitution 
coefficient a for the fluid-fluid particle collisions must approach unit as the mass ratio 
goes to zero. Hence, 

f* {2) (i>i,«2,qi2 = a) -»• xWVhcsMVhcsMi <Phcs(vi) = -^e - "*, (C.3) 

where is the equilibrium pair correlation function at contact. Then, the integrals 
in equation (jC.l|) can be performed, and the leading contribution to the cooling rate 
becomes 

4,r(d-i)/2 

c ''"~7ff(^ (1 - ah ° (F) - (c ' 4) 
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